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We consider a new type of magnetic tunnel junction, which con- 
sists of two ferromagnetic tunnel barriers acting as spin filters (SFs), 
separated by a nonmagnetic metal (NM) layer. Using the transfer ma- 
trix method and the free-electron approximation, the dependence of 
the tunnel magnetoresistance (TMR) on the thickness of the central 
NM layer, bias voltage and temperature in the double SF junction are 
studied theoretically. It is shown that the TMR and electron-spin po- 
larization in this structure can reach very large values under suitable 
conditions. The highest value of the TMR can reach 99%. By an ap- 
propriate choice of the thickness of the central NM layer, the degree 
of spin polarization in this structure will be higher than that of the 
single SF junctions. These results may be useful in designing future 
spin-polarized tunnelling devices. 

I. INTRODUCTION 

In the past few years, tunnel magnetoresistance (TMR) in magnetic tunnel junc- 
tions consisting of two ferromagnetic metal (FM) electrodes separated by a thin 
insulator (FM/I/FM) has attracted much attention due to its promising applica- 
tions in magnetic field sensors and magnetic random access memory [1,2]. The 
TMR in FM/I/FM junctions depends on the degrees of electron spin polarization 
of the two FM electrodes [3], but the lack of nearly perfectly spin-polarized current, 
and temperature dependence of the polarization, have limited the TMR in these 
junctions [4,5]. If we use half-metallic electrodes in the magnetic tunnel junctions, 
which are completely spin-polarized at T=0 K due to complications such as disorder 
and surface effects, half-metallicity is destroyed and thus the spin polarization and 
the TMR arc decreased [6]. However, by exploiting the spin filtering phenomenon, 
one can create 100% spin polarization and obtain a giant TMR. 

The spin filtering effect in a ferromagnetic semiconductor (FMS) has been demon- 
strated dramatically in spin-polarized tunnelling experiments [7], In these exper- 
iments, using an Al/EuS/Al tunnel junction, Moodera et al obtained 85% spin 
polarization for tunnel electrons. In another study with EuSe barrier junctions [8] 
and in the presence of an external magnetic field, they reported near 100% spin 
polarization. More recently, LeClair et al [9] from the combination of a spin filter 
(SF) barrier and a FM electrode, obtained a large TMR in an Al/EuS/Gd tunnel 
junction, which is a new method for spin injection into semiconductors [10]. When 
a FMS layer such as EuS is used as a tunnel barrier, due to the spin splitting of 
the conduction band in the FMS, tunnelling electrons see a spin-dependent barrier 
height. Thus the probability of tunnelling for one spin channel will be much larger 
than the other, and a highly spin-polarized current may result. Although the TMR 
using other SF barriers has also been studied [11-13], the results have shown only 
little success. 
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The TMR in single and double SF junctions has also been investigated theoret- 
ically in recent years. Based on the two band model and free-electron approxima- 
tion [14], Li et al [15] studied the tunnelling conductance and magnetoresistancc 
of FM/FMS/FM junctions. The results showed that a decrease or increase in tun- 
nelling current strongly depends on the magnetization orientations both in the elec- 
trodes and in the FMS layer. In another theory with double SF junctions, Worledge 
et al [16] studied the TMR of the NM (nonmagnetic metal) /FMS /FMS /NM junc- 
tions, using a simple model, and obtained a large magnetoresistancc In a recent 
paper, Wilczynski et al [17] investigated tunnelling in NM/FMS/NM/FMS/NM 
junctions in a sequential tunnelling regime. They found a strong enhancement of 
the TMR with increasing spin splitting of the barrier height in the ferromagnetic 
barriers. However, they have not investigated the effect of electric field and tem- 
perature on the tunnel currents and spin polarization. Thus, other aspects of this 
structure, such as the voltage dependence of TMR, remain to be explained. 

In the previous paper [18], using the single-site approximation for the 
NM/FMS/NM/FMS/NM double SF junction, we studied the effect of the thick- 
ness of the central layer on the spin polarization of tunnel electrons at T—Q K. We 
showed that the tunnelling spin polarization has an oscillatory behavior with the 
thickness, which is due to the spin asymmetry of the reflection at the FMS/NM 
interfaces. 

In this paper, using the transfer matrix method and the free-electron approxi- 
mation, we have extended our previous work [18] to investigate the effect of the 
thickness of the central layer, temperature and applied bias on the TMR for tun- 
nelling through a double SF junction. The paper is organized as follows. In section 
2, the transfer matrix approach of the spin-polarized tunnelling through the dou- 
ble SF junctions is described. In section 3, the numerical results for the TMR 
and the spin polarization of the tunnel currents are discussed. The conclusions are 
summarized in section 4. 



II. MODEL AND FORMALISM 



The system we consider here is a NM/FMS/NM/FMS/NM sandwich structure 
in the presence of an applied bias V a , which is depicted in Fig. 1. For simplicity, 
we assume the FMS layers, which act as SFs, are made of the same material, while 
the outer NM electrodes, which are semi-infinite, and the central layer are made of 
the same metal. This structure is grown in the x direction. In this case, in a free- 
electron approximation for the spin-polarized tunnelling electrons, the longitudinal 
part of the effective one-electron Hamiltonian can be written as 
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+ U j (x) + Vf , 
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where m* is the electron effective mass in the jth layer and 



Uj{x) 



0, 

-eV a x/(L - c) + U 2 , 
-eV a b/(L-c), 
-eV a (x - c)/(L - c) + U 4 , 
-eV a , 



x < , 
< x < b , 
b < x < b + c . 
b + c < x < L 
x > L , 



(2) 



where Ui and U4 are, respectively, the barrier heights of the left and right FMS layer 
at above Tc, b and d are the barrier widths, c is the width of the middle NM layer 
and L = b + c + d. The third term in Eq. (1) which is a spin-dependent potential, 
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denotes the s — f exchange coupling between the spin of tunnelling electrons and the 
localized / spins in the jth FMS layer. Within the mean field approximation, Vf is 
proportional to the thermal average of the / spins, (S z ) (a 7/2 Brillouin function), 
and can be written as V? = —Ia(S z ). Here, a = ±1, which corresponds to a =f , j, 
respectively and / is the s — f exchange constant in the FMS layers. 

The Schrodinger equation for a barrier layer under the influence of an applied 
bias can be simplified by a coordinate transformation whose solution is a linear 
combination of the Airy function Ai[p(z)] and its complement Bi[p(z)] [19]. Consid- 
ering all five regions of the NM/FMS/NM/FMS/NM junction shown in Fig. 1, the 
eigenf unctions of the Hamiltonian (1) with eigenvalue E x have the following forms: 



A 3a e ik3X + B 3a e~ ik3X , o 
A ia k\[p ia (x)] + B 4a Bi[p 4a (x)], b 



A 5 „e ik * x + Bj 



hot 



x < , 
< x < b , 
b < x < b + c . 
+ c < x < L 
x > L , 



(3) 



where the coefficients Aj a and Bj a are constants to be determined from the bound- 
ary conditions and 



h = y/2m\E x /h , 



(4) 



with 
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Although the transverse momentum ky does not appear in the above notations, 
the effects of the summation over ky will be considered in our calculations. 



A. Transmission coefficients 



By using the boundary condition such that the wavefunctions and their first 
derivatives are matched at each interface point Xj, i.e. tpj a {xj) = ipj + i, <T (xj) and 
(m*)~ 1 {d^j t(7 (xj)/dx} = (m* +1 )~ 1 {dtpj + i_ cr (xj)/dx], we obtain a matrix formula 
that connects the coefficients A\ a and Bi a with the coefficients A 5cr and B 5a as 
follows: 
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Since there is no reflection in region 5, the coefficient B^ a in Eq. (3) is zero. In 
this case the transmission coefficient of the spin a electron which is defined as the 
ratio of the transmitted flux to the incident flux, for the double SF structure shown 
in Fig. 1, can be written as 



Ta(E x ,V a ) = 



k^m\ 
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1 



M 11 
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where M^ tal is the left-upper element of the matrix M to tai which is defined in Eq. 
(11). Note that the transmission coefficient depends on the longitudinal energy E x , 
the applied bias V a and the spin orientation. 



B. Spin polarization and TMR 

The spin-dependent current density for single or double SF junctions at a given 
applied bias V a can be calculated within the free-electron model [20] : 

emlk B T f°° T( „ v]] \ l + cxp[(E F -E x )/k B T] \ 

Ta{E x , V a ) In < - — r-r= -Trrrr > dE x , (13) 



4rf Jo " ' 1 1 + exp[(£ F - E x - eV a ) /k B T] 

where k B is the Boltzmann constant, T is the temperature, and Ep is the Fermi 
energy. 

The degree of spin polarization for the tunnel current is defined by 

Jf + -h 

where Jj (Jj) is the spin-up (spin-down) current density. For the double SF junc- 
tion, this quantity can be obtained when the magnetizations of two FMS layers are 
in parallel alignment. 
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The tunnel conductance per unit area is given by G = J2a J<r/V a - In this case, 
the TMR can be described quantitatively by the relative conductance change as 

TMR= G II _G Ii; (15) 

where and G-[i correspond to the conductances in the parallel and antiparallcl 
alignments of the magnetizations in the FMS layers, respectively. 



III. NUMERICAL RESULTS AND DISCUSSION 



The numerical calculations have been performed for a NM/EuS/NM/EuS/NM 
double SF junction in which, for simplicity, we assume that the EuS layers have 
the same thickness b = d=0.5 nm. The appropriate parameters for EuS which have 
been used in this paper are: T c =16.5 K [21], 5=7/2, 7=0.1 eV [22], m* uS = 1.5 m e 
[23] and U2 = U4 = Ef+0.75 eV. In the NM layers, the electron effective mass and 
the Fermi energy are taken QbS 777 M — 777 q and Ep— 1.2b eV. Here m e is referred to 
the free-electron mass. In this study, we calculate the spin currents, tunnelling spin 
polarization and the TMR by using Eqs. (13)-(15), respectively. In our considered 
system, the magnetization orientation (i.e. the / spins' direction) in the left EuS 
layer stays fixed but the other EuS layer is free and may be switched back and 
forth by an external magnetic field (see Fig. 1). Thus, when the magnetizations 
of two FMS layers are parallel, spin-up and spin-down electrons see a symmetric 
structure, while for the antiparallcl alignment these electrons see an asymmetric 
structure. This structural asymmetry results from the difference in the two barrier 
heights for each spin channel. 

Figure 2 shows the TMR as a function of the thickness of the central NM layer at 
T=0, 0.45 and 0.9 Tq, when the bias voltage V a =50 mV is applied to the junction. 
It is obvious that the TMR oscillates with increasing thickness c and have well- 
defined peaks in which the TMR can reach 99% in some structures. The height of 
these peaks decreases with increasing c. The oscillatory behavior is related to the 
quantum-well states of the central NM layer and the spin-polarized resonant tun- 
nelling. On the other hand, due to the temperature dependence of spin splitting in 
the EuS layers, the barrier heights become spin-dependent, so that with decreasing 
temperature, this spin splitting, and thus the TMR, increases. 

To understand the physical origin of the TMR and the oscillations, we study the 
energy dependence of the transmission coefficients through the double SF junction. 
In Fig. 3, we have plotted the spin-dependent transmission coefficients at T—0 
K for c=0.50 nm, which corresponds to a flat area between the peaks, and for 
c=0.72 nm, which corresponds to a local maximum in the TMR. Because of the 
quasibound states in the central NM layer, the transmission coefficients reach unity 
at the resonance peaks which become sharper in the low incident energy region, 
since in this energy region the resonance levels are more strongly quantized. The 
results for c=0.50 nm show that there is one resonance level in the quantum well 
for both spin orientations and magnetic alignments. All these resonance levels 
are far from the Fermi energy. However, the transmission coefficient for spin-up 
electrons in the parallel alignment is higher than the antiparallel alignment, and 
for spin-down electrons in both alignments. Thus, there is a small difference in 
the current density J (=J| + Jj.) in both magnetic configurations, which gives 
rise to relatively small TMR at this thickness of the central layer. For c=0.72 
nm the resonance states shift to the lower energy side and a new resonance level 
slightly below Ep appears, which is active only for spin-up electrons in the parallel 
alignment. Therefore, there is a large difference in the current density J in both 
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alignments and consequently gives rise to large TMR. It is clear from the Fig. 
3(b) and 3(d) that, for each thickness c, the resonance levels for spin-up and spin- 
down electrons in the antiparallcl alignment coincide. The reason is that, for the 
antiparallel alignment, electrons with up (down) spin are easy (difficult) to tunnel 
into the central NM layer, and difficult (easy) to tunnel out of it; thus, both the 
spin-up and spin-down electrons see the same effective height of the barriers during 
the tunnelling process through the whole system. 

In Fig. 4 we show the TMR as a function of the applied bias at T=0 K for the 
thicknesses c correspond to Fig. 3. With increasing bias voltage, the TMR for 
c=0.50 nm decreases very slowly because, in this range of the applied voltage, the 
discrepancy between the conductance for the parallel alignment and that for the 
antiparallel alignment increases only slightly. On the other hand, the results show 
that, for c=0.72 nm at the beginning, the TMR slowly decreases with increasing the 
bias voltage. However, for the voltages higher than V a =80 mV, it quickly decreases. 
The reason is that, at higher voltages, one of the resonance levels becomes active 
for both magnetic alignments which drastically reduces the TMR. 

It should be noted that, for very low values of the applied bias and the incident 
energy, a numerical instability is occurred in some of our calculations, which is 
due to the use of exact Airy functions. This instability is overcome by using the 
asymptotic forms of Airy functions [19] and numerical analytical methods. 

In Fig. 5 we have shown the spin polarization of the tunnelling current versus 
normalized temperature T/Tc for single and double SF junctions to reveal the SF 
effects of the FMS layers from another point of view. At high temperatures T > Tc, 
there is no spin splitting E ex in the conduction band of the EuS layers and the trans- 
mission coefficients for two spin channels coincide. Thus there is no TMR and spin 
polarization effect. As the temperature decreases from Tc, the barrier heights for 
spin-up electrons are lowered, while they are raised for spin-down electrons. This 
temperature dependence of the barrier height, which is attributed to the exchange 
splitting of the EuS conduction band, greatly increases the tunnelling probabil- 
ity for one spin channel and reduces it for the other. In the parallel alignment, 
the tunnel current for spin-up electrons is much higher than the spin-down ones, 
which gives rise to the TMR and spin polarization effect. On further decreasing 
the temperature, this spin splitting and hence the difference in the barrier heights 
increases. Therefore, the TMR and spin polarization reach the highest values at 
T=0 K. The highest value of the spin polarization for the single SF junction can 
reach 77%, which is qualitatively in agreement with the experimental measurements 
[7] and the theoretical results [24,25], while for the double SF junction (in the par- 
allel alignment), it approaches 66% for c=0.50 nm and 99% for c=0.72 nm [18], 
which is due to the change in the positions of the spin-polarized quasibound states 
in the quantum well. Therefore, one can see that, for the double SF junctions, the 
TMR and the spin polarization of the tunnelling current can be controlled by the 
thickness of the central NM layer, the temperature, and the applied bias. 

In this study, the numerical results obtained for the TMR and the degree of 
spin polarization can be compared with the result of Wilczynski et al [17]. As we 
discussed above, the enhancement of spin splitting E ex in the conduction band of 
the FMS layers can be obtained, when the temperature decreases. In this case, the 
TMR and the spin polarization increases, and this behaviour is in agreement with 
the result obtained by Wilczynski et al in the sequential tunnelling regime. 
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IV. SUMMARY 



Based on the free-electron model, we presented a transfer matrix method for spin- 
polarized tunnelling through the double SF junctions. The effect of the quantum 
size, applied bias and temperature on the spin filtering and the spin transport 
process in the FMS/NM double junctions are examined. Numerical results indicate 
that the TMR oscillates as the thickness of the central NM layer increases. It 
is further confirmed that, for some thicknesses of the central NM layer, the spin- 
polarized resonant tunnelling can gives rise to large values for the spin polarization 
and the TMR, even at high temperatures (T < Tc). Therefore, in the system 
presented here, it is able to select an appropriate applied bias, temperature and the 
thickness to achieve a maximum TMR and spin polarization. 

Although the temperature for observing a 99% TMR is very low and the findings 
of the paper are not directly applicable to spintronics technology, the results may 
be useful in designing future spin-polarized tunnelling devices [26] . 

In the present model, we have neglected the generally important complications, 
such as the interface roughness, electron-electron interaction, magnetic-domain 
wells, / spin correlation, etc. These effects can play important roles in the spin 
transport process and reduce the efficiency of spin filtering. 
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Figure 1 



A. Saffarzadeh 



FIG. 1. Spin-dependent potential profile for NM/FMS/NM/FMS/NM tunnel junctions 
in the presence of a positive bias V a - The broken line in the FMS layers represents the 
bottom of the conduction band at T > Tc- Below Tc, due to the exchange splitting E ex 
in the conduction band, the barrier heights become spin-dependent, as indicated by the 
thin arrows for spin-up and spin-down electrons. The direction of magnetization in the left 
FMS layer is fixed in the +z direction, while the magnetization in the right FMS layer is 
free to be flipped into either the +z or — z direction, as indicated by the filled and hollow 
arrows, respectively. 
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Figure 2. 



A. Saffarzadeh 



FIG. 2. TMR as a function of the thickness c in the EuS/NM double SF junction at 
T=0, 0.45 and 0.9 T c . 
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A. Saffarzadeh 



FIG. 3. Spin-dependent transmission coefficients, In Tc(E x ), as a function of energy E x 
the EuS/NM double SF junction for c=0.50 and 0.72 nm at T=0 K. 
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FIG. 4. TMR as a function of applied bias V a in the EuS/NM double SF junction for 
c=0.50 and 0.72 nm at T=0 K. 
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FIG. 5. Electron-spin polarization P as a function of normalized temperature T/Tc for 
the EuS/NM single and double SF junctions. 
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